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ABSTRACT 

We study the cross-correlation of distribution of galaxies, the Sunyaev-Zel’dovich (SZ) 
and X-ray power spectra of galaxies from current and upcoming surveys and show 
these to be excellent probes of the nature, i.e. extent, evolution and energetics, of 
the circumgalactic medium (CGM). The SZ-galaxy cross-power spectrum, especially 
at large multipoles, depends on the steepness of the pressure profile of the CGM. 
This property of the SZ signal can, thus, be used to constrain the pressure profile of 
the CGM. The X-ray cross power spectrum also has a similar shape. However, it is 
much more sensitive to the underlying density profile. We forecast the detectability 
of the cross-correlated galaxy distribution, SZ and X-ray signals by combining South 
Pole Telescope-Dark Energy Survey (SPT-DES) and eROSITA-DES/eROSITA-LSST 
(extended ROentgen Survey with an Imaging Telescope Array-Large Synoptic Survey 
Telescope) surveys, respectively. We find that, for the SPT-DES survey, the signal- 
to-noise ratio (SNR) peaks at high mass and redshift with SNR ~ 9 around Mh ~ 
10 13 /i _ 1 Mq and 2 ~ 1.5-2 for flat density and temperature profiles. The SNR peaks at 
~ 6(12) for the eROSITA-DES (eROSITA-LSST) surveys. We also perform a Fisher 
matrix analysis to find the constraint on the gas fraction in the CGM in the presence or 
absence of an unknown redshift evolution of the gas fraction. Finally, we demonstrate 
that the cross-correlated SZ-galaxy and X-ray-galaxy power spectrum can be used as 
powerful probes of the CGM energetics and potentially discriminate between different 
feedback models recently proposed in the literature; for example, one can distinguish 
a ‘no active galactic nuclei feedback’ scenario from a CGM energized by ‘fixed-velocity 
hot winds’ at greater than 3er. 
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1 INTRODUCTION 

Galaxy formation theories predict that galaxies form by the 
collapse of dark matter haloes when the density crosses a 
critical overdensity. Baryonic matter falls into the gravita¬ 
tional potential driven by dark matter, cools down, forms 
stars and hence galaxies. The standard scenario predicts 
that the amount of baryons in a galactic halo should approxi¬ 
mately be a constant fraction (~ 16%) of the total halo mass 
(dark matter + baryonic matter). This baryon-to-total halo 
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mass fraction is known as the cosmic baryon fraction (CBF). 
On the other hand, observations have detected only a small 
fraction of this CBF. More than half of the baryons of CBF 
are missing f rom the galaxies according to soft X-ray ab¬ 
sorption line dAn derson fe Breg manll2010l ) and emission line 
searches dMiller fc Bregman 20151 ) in galactic haloes. This is 
an important problem for all galaxy formation studies, for 
with more than 50% of the baryon reservoir unaccounted for, 
the sources and sinks for forming stars become correspond¬ 
ingly uncertain. Recent detection of substantial amounts of 
cooler gas within the halo virial radius f ails to significantly 
alleviate this problem (IWerk et al.ll2014l ). 
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Although numerical simulations suggest the presence of 
hot coronal gas, it has been difficult to detect this gas due 
to the faintness of its X-ray emission. However, there are 
indeed recent observations which indicate the presence of 
significant amounts of hot coronal ga s extended over a large 
regio n around massive spiral galaxies (lAnderson fc Bregman 
2011); Dai et all 120121 ; lAnderson et al.l 20131 : Bogdan et al. 
2C)13a|07" Also, the inferred ram pressure stripping of 


satellite galaxies around the Milky Way Galaxy supports 
the presence of the hot halo gas with nearly flat den- 
sity (n ~ 10~ 3,5 cm~ 3 ) out to large galacto-centric radius 

([Grcevich fe Putmanll2009l: IPutman et al.ll2012l: iGatto et al.l 

l2013lh This gas, known as the circumgalactic medium 
(CGM), may account for some of the missing baryons from 
the galaxies. 

The CGM is the gas surrounding the central, optically 
visible part of the galaxy within its host dark matter halo. 
This is the bridging medium that connects the interstellar 
medium (ISM) to the intergalactic medium (IGM). During 
galactic evolution, the galaxy accretes matter from its sur¬ 
rounding IGM and also ejects some material in the form of 
galactic winds resulting from feedback processes like super¬ 
novae (SNe) and active galactic nuclei (AGN). The CGM, 
being the intermediate medium, is most affected by these 
processes and may contain important clues about galaxy 
evolution, making it a promising tool for the study of the 
processes affecting the galaxy evolution. 

The presence of the hot halo gas surrounding the galax¬ 
ies, in the form of the CGM, causes the Sunyaev-Zel’dovich 
(SZ) distortion of cosmic microwave background radiation 
(CMBR). The thermal SZ (tSZ) distortion of the CMBR 
due to the CGM is small comp ared to the tSZ distor¬ 
tion caused by the galax y clusters (IPlanck Collaboration Xlj 
120131 : Isingh et aPl2015l l. The detectability of the tSZ signal 
from a system can be enhanced by cross-correlating this sig¬ 
nal with another signal ori ginati n g from the same source like 


the distributio n of haloes dF ang et al.||2012f) or the gravita - 


tional lensing (IVan Waerbeke et al.l 20141 : ilia et al] 


We cross-correlate the tSZ signal from the CGM with the 
distribution of the galaxies. The cross-power spectrum be¬ 
tween the SZ signal and distribution of galaxies can be 
thought of as the SZ-galaxy cross power spectrum. It can 
be obtained by combining a high resolution a CMB survey 
such as SPT (South Pole Telescope) with an overlapping 
optical survey. 

In addition to the SZ-effect, the hot CGM also manifests 
itself in X-ray emission through bremsstrahlung. Combining 
the X-ray observations with optically selected galaxies can 
give the X-ray-galaxy cross power spectrum, which enhances 
the detectability of the X-ray emitting gas. The X-ray emis¬ 
sion from the CGM is more sensitive to the underlying gas 
distribution than the SZ effect and it also breaks the de¬ 
generacy between the gas density and temperature which is 
present in the SZ effect. Here, we study the prospects for the 
cross-correlation of the soft X-ray emission from the CGM 
with the distribution of galaxies. This can be used as an 
additional probe to constrain the properties of the CGM. 

We also compute the X-ray-SZ cross-power spec¬ 
trum for the CGM. Significant effort has been made 
to forecast/detect the X-ray-SZ cross-correlation on large 
scales by cross-correlating the CMB maps generated by 
WMAP (Wilkinson Microwave Anisotropy Prob e)/Planck 


surv eys with the ROSAT (Rontgen Satellite) all sky sur¬ 
vey (iDiego et al' 20031: Hernandez-Monteaeudo et al.ll2004l . 


I2OO6I : iHinshaw et al. 120071 : lHaiian et al.ll2013h . The X-rav- 

SZ cross-correlation measured for the galaxy clusters is par¬ 
ticularly useful to constrain the cosmological parameters as 
the number of clusters strong l y dep ends on the underlying 
cosmology dHurier et al.ll2014l . 120151 1. We look into the pos¬ 
sibility of detecting the contribution of the CGM to the X- 
ray-SZ cross-power spectrum with SPT and eROSITA (ex¬ 
tended ROentgen Survey with an Imaging Telescope Array) 
combination. 

Simulations suggest that feedback processes are re¬ 
quired to avoid the over-cooling problem and formation of 
excessively massive galaxies (in terms of stellar mass). A 
number of feedback mechanisms have been proposed that 
reproduce many observed galaxy properties despite having 
different implementations and physical motivations behind 
them. The CGM, thus can provide additional constraints on 
these simulations as the CGM properties are la rgely affected 
by the variation in the feedback mechanism dSuresh et al.l 
l2015ll 

This paper is organized as follows: In section-0 we de¬ 
scribe the SZ-effect. In section-[3] we estimate the cross¬ 
correlation between the SZ-effect from the hot CGM and 
the distribution of galaxies and forecast the detectability 
of the SZ-galaxy cross power spectrum. In section-[f] we de¬ 
scribe the X-ray emission from the hot CGM and forecast 
the detectability of X-ray-galaxy cross power spectrum. In 
section-[5] we compute the X-ray-SZ cross power spectrum. 
In section-|6] w e forecast the constraints on the CGM proper¬ 
ties. In section-|T]we discuss the possibility of differentiating 
between various feedback models. In section-[8] we conclude 
by summarizing our main results. 


2 SZ DISTORTION FROM HOT GALACTIC 

HALO GAS 

For simplicity, we assume that the mass fraction of the 
CGM is independent of the host halo mass. Observations 
indicate that the fractional mass in the s tellar compo- 
nent of galaxies is ~ 0.05 llMo et al.l Il998l: iDutton et al.l 

l20ld : iLeauthaud et al1l20ld : iMoster et alJ2Qloh . This corre- 

sponds to a fractional mass in gas, / gas = 0.11, as the total 
fractional mass in baryons in a galaxy is ~ 0.16. Due to the 
uncertainty in the amount of the CGM, we also calculate 
some key results with a smaller gas fraction, / gas =0.05. We 
assume that the gas is uniformly distributed in the galactic 
halo with a temperature given by the virial temperature of 
the halo. We also show the effect of different density profiles 
of the CGM on its cross power spectrum. 

The cosmological parameters that we have used are 
driven by the joint analysis of CMB anisotropies along with 
observations from the Baryon Acoustic Oscillations surveys. 
The resulting best fit cosmological parameters, especially 
ag, are in tension with those obtained from galaxy clus¬ 
ter and weak lensing studies. In particular, one can have a 
mismatch in the galaxy cluster counts by a factor of two 
due to the difference in the adopted value of as- Although, 
much efforts have been ma de to calibrate cluster masses 
(IPlanck results XXIVl l2015l l. crucial for doing cosmology 
with clusters, the simplest explanation for the mismatch can 
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be a remaining systematic mass bias at the galaxy cluster 
scales. Similarly, unknown systematics can also lead to the 
amplitude of the fluctuation spectrum inferred from weak 
gravitational lensing to be lower than that inferred from 
CMB. However, this tension has been lifted to a certain ex¬ 
tent by the ‘first’ cosmological results from the Dark Energy 
Survey (DES, which is one of the surveys that we consider 
in our analysis), where th e estimated o ~8 is consiste nt with 
the Planck measurement (IDES collaboration! l2015l l. There 
are recent indications that masses estimated from velocity 
dispersions (for galaxy clus ters) may have m ore robustness 
than previously envisaged dRines et al.ll2015h . Taking posi¬ 
tively these developments, which bring increased consensus 
among different cosmological r esults, we adopt the Pla nck 
CMB cosmological parameters (IPlanck results XIIll2015h as 
our fiducial choices. We do comment on the impact of cos¬ 
mological parameters/degeneracies on our results in section- 

EH1 


2.1 Thermal SZ effect 


The inverse Compton scattering of the CMB photons by 
high energy electrons distorts the C MB spectrum giving rise 
to the tSZ (ISunvaev fe Zebdovich Ill969l f. The tSZ effect is 
represented in terms of a dimensionless parameter, known 
as the Compton y-parameter, defined as y = f dl ttTeT * e i <rT f 
which for a flat density and temperature profile becomes 
y = ( kbT e n e aTL)/(m e c 2 ) where or is the Thomson scat¬ 
tering cross section, T e is the gas temperature ( T e » T 7 ), 
n e is the electron density of the scattering medium, and L 
is the distance travelled by the photons through the scat¬ 
tering medium. The electron density n e = P9as of the gas 
is determined by the condition that the total hot gas mass 
within the virial radius is a fraction f g = 0.11 of the total 
halo mass. For a halo of mass M at redshift z, the virial 
radius is given by 


R vir = 0.784 h ~ 3 



M 

10 8 h" 1 



O.M A (z) 
Q.m{z) 18 - 7 T 2 



where A(z) = 187 r 2 + 82 d — 39 d 2 is the critical overdensity 
with d = SIm(z) — 1 and Qm{z) = Qm (1 + z) 3 /E 2 (z), where 
E{z) = + Om( 1 + z ) 3 ■ 


The angular power spe ctrum of the tSZ effect 
(Komatsu fo Kitavamall 19991) is given by 


(JVV — (j yv ' 


+ cr- 


(4) 


where Cf v,lh is the 1-halo or Poisson term and C/ y,2h is the 
2-halo or clustering term. These two terms can be written 


as 
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where r{z) = (1 + z)Da is the comoving distance, is 
the differential comoving volume per steradian, P m {k,z) is 
the matter power spectrum, b{M, z) is the linear bias fac¬ 
tor dSheth fe Tormenlll999lf . dn ^’ z - ) is the differential mass 
function, g{x v ) = x„ coth(a:„/2) — 4 is the frequency de¬ 
pendence of the tSZ effect with x v = k t C j MB ■ We compute 
all the power spectra in dimensionless units throughout the 
paper. All calculations in this work are done in the Rayleigh- 
Jeans (RJ) limit {g{x v ) —> —2). The term W y {z) present in 
the 2 -halo term is defined as 

r M max 1 

W y {z) = j^ dM—(M,Z)b(M,z) yi (M,z) ( 6 ) 


We use the Sheth-Tormen mass function given by 


dn 
dM " 


dM 


/2aV 2 p m _- a n 2 /2\ 

r d log (t I 

1 7T M 2 ® 

L d log M \ 


x j^l + ^ at j dM , 


(7) 


where A = 0.322184, a = 0.707, p = 0.3 dSheth fe Tormenl 
l 200 lh and v — D where 5 C = 1.68 is the criti¬ 

cal overdensity, D g {z) is the growth factor and a(M) is 
present day smoothed (with top hat filter) variance. We 
take Zmax = 8 as the upper redshift integration limit and 
Mmax = 10 13 /i _1 Mg as the upper mass integration limit. 
The lower mass integration limit is set by the condition that 
the gas cooling time-scale is large r than the halo des truction 
time-scale (explained in detail in lSingh et al.l (120151 ')'). 


3 SZ-GALAXY CROSS-CORRELATION 
3.1 The angular power spectrum for tSZ 


The angular Fourier transform of the Compton y-parameter 
(under flat sky approximation) is 


Vi 


47 tR s 


f 

Jo 


dxx 2 y3D(x) 


sin(lx/l s ) 
{lx/l s ) 


( 2 ) 


where x = r/R s is the dimensionless scaled radius, R s = 
R v ir/c(M, z) is the scale radius , l s = dA(z) IRs and c(M, z) 
is the concentration parameter (iDuffv et al.ll2008h . To calcu¬ 
late the SZ-effect from the CGM, we truncate the integration 
in equation-® at r = R v i r ■ The 3D-radial profile 2/3 d{x) is 


= - -n e (x)kbT e ( x) (3) 

m e c 


3.2 The angular power spectrum for the 
distribution of galaxies 

For simplicity, if we assume that the mass and red¬ 
shift of a galaxy can be measured accurately, the prob¬ 
ability that a given galaxy lies in the a th redshift 
bin Zobst [Zobs,min, Zo bs ,max] and b th mass bin M oba t 
[Mobs,min j M^bs mnx] is represented by a selection function 
defined as (lOguri fc Takadall201ll : iFang et aPl2012l ) 

Sab{M, Z) Q(z Z 0 b s ,ruin')®{Zobs,max z) 

X 0(M - M^bs,min)®{Mobs,max ~ M) (8) 

where 0 is the Heaviside step function. The power spectrum 
for the distribution of galaxies/galactic haloes in the {ab) th 
and {a'b') th bins is 

c% ah , a , v) = f dz^Pm(k = J-)wZ b {z)wZb'{z) (9) 


2/3 d{x) 
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Figure 1. The cross power spectra of the SZ effect in the RJ limit and the distribution of galaxies in different mass and redshift bins. Here 
green (dashed) lines, blue (dot-dashed) lines and red (solid) lines represent the 1-halo term, 2-halo term and the total signal respectively. 


with Wab(z) defined as 

wUz) = -^of dM-^S(M,Z)S ab (M,Z)b(M,Z) (10) 

where is the two-dimensional angular number density 
of the galaxies in ( ab) th bin and is given by 

nl? = I dz ^ J dMS ab (M,Z)-^(M,Z) ( 11 ) 


3.3 SZ-galaxy cross-correlation power spectrum 


The SZ-galaxy cross power spectrum is the cross-correlation 
between the SZ signal and the distribution of galaxies. For 
the galaxies in the ( ab) th bin, the SZ cross power spectrum 
is 


s~iyh _ s-iyh,lh . s^yh,2h 

°Z,(a&) — °Z,(ab) ^i^ a b) 


( 12 ) 


where and CfJ 1 ’?! 1 are the 1-halo and 2-halo 

l,{ab) l,{ab) 

terms respectively and these terms c an be written as 
llOguri fe Takadall201ll : iFang et ahll2012l ') 


y-iyh,lh _ 

°Z,(ab) — 


(13) 


l d ^P™(k = ^) W Uz) W ?(.z) ( 14 ) 

Due to the presence of S ab (M, Z) in Ea 1 131 and 1 141 only the 
galaxies lying in the a th redshift bin contribute to the SZ 


cross power spectrum. For the 2-halo term, even the galaxies 
lying outside the b th mass bin contribute because of the the 
presence of W?(z) in Eq[Q](see Eq(S]). 

In Figure- 0 we show the SZ-galaxy cross power spec¬ 
trum for a few mass and redshift bin combinations. The 
galaxies are binned in mass (total halo mass) and redshift 
with A log(M) = 0.2 and Az = 0.1 respectively. The mean 
redshift increases from left panel to right panel and the mean 
halo mass increases from top panel to bottom panel. Most 
of the contribution to the total cross power spectrum comes 
from the 1-halo term. 


3.4 CGM density profile and the cross-power 
spectrum 

For a flat pressure profile for the CGM, the SZ-galaxy cross 
power spectrum shows oscillations if l > Imax , where l ma x 
depends on the mass and redshift of the galaxy. In Figure- m 
we show lmax as a function of the mean halo mass for differ¬ 
ent redshift bins. The oscillations begin when the multipole 
l corresponds to an angular size ~ |x virial radius of the 
galaxy. The reason for these oscillations is the truncation of 
the signal at R v i r - 

The shape of the SZ-galaxy cross power spectrum is 
sensitive to the pressure profile of the CGM. Since we have 
fixed the temperature of the CGM, the density profile of 
the CGM can be constrained using the cross-correlated 
SZ-galaxy power spectrum. In order to estimate the effect 
of different density profiles on the SZ cross power spectrum, 











CGM Cross-correlations 5 


we parameterize the density profile by 7 gas , defined by 
pgas oc [l + (r/i? s ) 7gas ]' _1 In Figure-Q we show the total SZ 
cross-power spectrum for mass bins:[10 12 ' 8 , 10 13 '°](/i - 1 M©) 
and redshift bins:[0.4, 0.5] for the following three density 
profiles, 


Profile-(a): 7 gas =0 =>• a flat density profile. 
Profile- (b): 7 gas =1 =>■ p gas oc [1 + (r/f? s )] _1 
Profile-(c): 7 gaa =3 => p gas oc [1 + (r/7? s ) 3 ] -1 


From profiles-(a) to (c), the density becomes steeper 
and more centrally concentrated. Above density profiles 
are similar to the /3-model (with /3 = 2/3) with a central 
core followed by a gradual decrease in the density. This 
choice of the density profile is inspired by the observation 
of nearly flat di stribution of the hot halo g a s in the M ilky 
Way Galaxy dGrcevich fe Putmanl I200S : IPutman et al.l 
20j^ jGattoet id! _ 2013l l an d the simul at ion r esults of 


Le Brun et al. ( 2015l l and ISuresh et al.l (l2015l l , which 


predict a flatter distribution of the gas compared to the 
Navarro-Frenk-White (NFW) profile in galaxies due to the 
presence of feedback processes. In all 3 cases the SZ signal 
is truncated at the virial radius. As the density profile 
becomes steeper from (a) to (c) the value of Imax shifts 
from 12000 to 14000 whereas there are no oscillations in 
the case of profile-(c). This shift occurs since with the 
steepening of the density profile, the pressure at the virial 
radius decreases. 

In Figure-® at small /-values (~ 3000), the SZ cross 
power spectrum for the three density profiles are almost 
identical but there is a significant difference between the 
profiles near / ~ 10 4 . This is because the steepening of the 
density profile increases the power at small angular scales or 
large / values. The value of / where the shape of the cross¬ 
power spectrum is significantly different for different density 
profiles depends on the mean mass and redshift of the bin 
in a similar way as for l m ax■ Therefore, the shape of the SZ- 
galaxy cross-power spectrum at these /-values can be used 
to determine the slope of the density profile of the CGM. 
However, the use of this method is limited by the resolution 
of the CMB survey. Presently the SPT survey has the best 
resolution and it goes up to / ~ 10 4 . High mass and low 
redshift galaxies are better choices for this purpose as the 
SZ power spectrum for different profiles is distinguishable at 
l < 10 4 for these galaxies. 


3.5 Detectability of the SZ-galaxy 
cross-correlation signal 

3.5.1 Cumulative signal-to-noise ratio 

Given a survey one can estimate the detectability of the SZ- 
galaxy cross-correlation signal for galaxies binned according 
to their masses and redshift using the cumulative signal-to- 
noise ratio (SNR) defined as 

(|) 2 = (15) 

where M^ {ab) is the covariance matrix which represents the 
uncertainty in the measurement of Cf^ ab y The covariance 



Figure 2. Imax is shown as a function of the mean halo mass 
for redshift bins 2 : : [0.4-0.5] (red solid line) , 2 : : [0.7-0.8] (green 
dashed line), 2 : : [1.0-1.1] (blue dotted line) and 2 : : [1.3-1.4] 
(brown dot-dashed line). 



/ 


Figure 3. The cross power spectra of the SZ effect (in 
RJ limit) and the distribution of the galaxies for M ^ : 

[10 128 ,10 13 , 0 ](/i _ 1 Mq) and 2 : [0.4, 0.5] for the density profile- 
fa) (red solid line), profile-(b) (green dashed line) and profile-(c) 
(blue dot-dashed line) i.e. 7 ga s = 0, 1 & 3 respectively. In the in¬ 
set, we show the corresponding density profiles. 


matrix is given by 


Mi 


yh 


ll',(ab,a'b') 


$11' 

fsky{2l + 1)AZ X 


\AyyAhh 
ly-'Z ^ l,(ab,a 



(16) 

where f a k y is the fractional sky coverage for a survey, A l is 
the 1-bin size used to calculate the power spectrum and Ci’s 
represent the power spectrum including the noise contribu¬ 
tion (i.e. Ci + N }) where i stands for yy, hh and yh. For 
simplicity, we neglect the non-Gaussian contribution for the 
/-range considered: this is particularly true for the massive 
galaxies. Note that the dominant contribution to the noise 
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Table 1. Specifications of surveys 


Survey 

n s 

(deg 2 ) 

J sky 

V 

(GHz) 

#FWHM 
(■ arcmin ) 

(7 t 

(l*K) 




95 

1.6 

26.3 

SPT 

2500 

6% 

150 

l.l 

16.4 




220 

1.0 

85 


Survey Q s 

(deg 2 ) 


Jsky 


z range 


DES 

5000 

12% 

0.1-1.4 

LSST 

20000 

48% 

0.1-1.4 


Survey 

(deg 2 ) 

J sky 

#FWHM 

(arcsec) 

Exposure time 
(ks) 

eROSITA 

All sky 

100% 

20 

2 


in the SZ-galaxy as well as X-ray-galaxy cross-power spec¬ 
trum comes from the instrumental noise. Thus, the inclusion 
of clusters has negligible effect on the error calculations and 
hence the SNR. Therefore, we neglect this contribution in 
our calculation. 

The instrumental noise simply adds to the SZ power 
spectrum and for a given CMB survey, it is given by 


N, 


yy,(CMB) _ 
l 


E kWksiB?, 


(17) 


where the summation in Eq |17l is over the different frequency 
channels, w = (ctt#fwhm/Tcmb) -2 where or is the rms 
instrumental noise per pixel, #fwhm is the full width half¬ 


maximum of the beam, Bf = exp|^—Z(Z + l)0| WHM /(8Zn2)j 


is the fourier transform of beam profile and s = —g{x v )/2 
is to rescale the result in RJ limit. The shot noise in the 
galaxy distribution power spectrum is 


N, 


hh,(g) 

( ab,a'b') 


2D tiaa'^bb' 

,L ab 


(18) 


For the cross-power spectrum NJ h = 0 i.e. Cf h = Cf h 
as the distribution of galaxies is not correlated with the in¬ 
strumental noise in the CMB surveys. 


3.5.2 CMB survey 

To detect the SZ-galaxy cross-correlation signal, we need 
a galaxy survey and a CMB survey with overlapping sky 
coverages. Since the SZ signal from the CGM becomes non- 
negligible compared to other contribution to CMB distortion 
at large Z-values ( at Z > 3000 ), we consider the SPT survey 
for this work. Presently, the optical survey which overlaps 
with the SPT survey, is the DES. The specifications of these 
surveys are given in Table- ©• 

Note that, the dominant contribution to the covariance 
in EqJT6] comes from Nj JV ’ < ' CMB ' > /n 2 a B for Z > 3000. The 
contribution from other terms ( C^ v C^ h , N^ y:( - CMB \ 

Cr/n% and C? h C? h ) is negligible compared to this 
term. For example, for the haloes in the mass bin 
Mh : [10 12 ' 8 , 10 13 '°](/i _1 Mq) and redshift bin a : [0.4, 0.5], 



Figure 4. Contours of the cumulative signal-to-noise (SNR) for 
the measurement of the cross-power spectrum of the SZ effect 
with the distribution of galaxies for the combined SPT-DES like 
surveys. The red (solid), green (dashed), blue (dotted) and brown 
(dot-dashed) lines represent the SNR 9, 6.5, 5 and 3 respectively. 
The upper redshift limit plotted here is more than the highest 
redshift probed by DES. 


the instrumental noise is nearly three to four orders of 
magnitude larger than Cf h . 

In figure-© we show the contours of cumulative SNR of 
the SZ-galaxy cross power spectrum of the CGM for the flat 
density profile. The galaxies are binned in mass and redshift 
with A log (M) = 0.2 and A z = 0.1 respectively. Here we 
take Imin = 3000. Note that the DES survey goes only up 
to t ~ 1.4 whereas in Figure- © we have shown the results 
up to z ~ 3 to show the decreasing contribution from high 
redshift haloes. The SZ signal increases with increasing halo 
mass for a given redshift as the amount of gas causing the 
SZ-effect increases with increasing halo mass. This results 
in higher SNR for higher mass galaxies compared to the low 
mass haloes in the same redshift bin. For a given halo mass, 
the SNR first increases with increasing redshift, reaches a 
maximum value and then decreases. As a combined effect 
there is an optimum spot with SNR ~ 9 at the high mass 
end of galaxies, at redshift around z ~ 1.5-2. Even for low 
mass galaxies (Mh ~ 2-4 x 10 12 /i _1 Mq) the SNR is ~ 3 in 
the redshift range covered by the DES survey. This makes 
the SPT-DES survey a useful tool to study the CGM in the 
low mass galaxies which is otherwise difficult to detect at 
higher redshifts. 

In the calculation of the SNR in this section as well 
as the constraints on the model parameters from the SZ- 
galaxy cross-power spectrum in section- 16711 andl7l we neglect 
the contamination from other astrophysical sources (kinetic 
SZ (kSZ), Cosmic infrared background (CIB), point sources 
etc.). However, the tSZ signal has a distinct frequency de¬ 
pendence and the relative contribution from other sources 
can be minimized using multi-frequency data from surveys 
like SPT. 
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Figure 5. The cross power spectra of the X-ray emission from the CGM and the distribution of galaxies in different mass and redshift 
bins. Here green (dashed) lines, blue (dot-dashed) lines and red (solid) lines represent the 1-halo term, 2-halo term and the total signal 
respectively. 


4 X-RAY-GALAXY CROSS-CORRELATION 


The hot CGM causing the SZ-effect also emits in the soft 
X-rays. This X-ray emission from the galaxies can be cross- 
correlated with the galaxy distribution to increase its de¬ 
tectability. The X-ray surface bright ness due to the prese nce 
of the hot gas in a direction 9 (eg. ICheng et alj (120031 )) is 
given by 

Mi W / w ' A(T ’ z)d * (19) 

where A(T, Z) is the cooling function which depends on the 
metallicity and temperature of the gas. We assume that the 
halo gas has a metallicity ~ 0.1 Z p > and use the cooling 
function from ISutherland fe Dopital (1 1993T) . For the metal¬ 
licity and mass range of interest, we calculate the emission 
from the CGM in the soft X-ray band (0.5-2.0 keV) assum¬ 
ing that the fraction of total energy in soft X-ray band is 

~ (exp 7 =t T —exp *iT r; ) , where Ei is the lower and E 2 
is the higher energy limit of the soft X-ray band. The mass 
and redshift range considered here corresponds to the tem¬ 
perature range ~ 10 6 -10' K. Therefore, these galaxies lie 
near the lower energy limit of the soft X-ray band used for 
this study. 

The fluctuations in the X-ray background in the direc¬ 
tion 9 are 


s{») 


S(0) 

< 5’sxrb > 


( 20 ) 


where < Ssxrb > is the mean surface brightness of the 


soft X-ray backgroun d (SXR B) and it inc ludes all possible 
sources of SXRB (see I Merloni et akl J2012I )). 


4.1 The X-ray angular power spectrum 

In analogy with the SZ-effect, the angular Fourier transform 
of fluctuations in the SXRB (for Z > 0) is given by 


si 


4-kRs 

IT 


poo 

/ dx X 2 S3d(x) 

Jo 


smilx/ls 

(lx /Is) 


where the 3-D-radial profile S3d(x) is 

1 n 2 (x)A(T, Z) 


S3D{x) = 


47t(1 + z) 4 < S SXRB > 


( 21 ) 


( 22 ) 


The X-ray power spectrum can be obtained by replacing 
g{x v )yi by si in Eq-©. 


s-jxx,lh _ 

r* 

dz^ 

v 't 

Jo 

dz 

s^yXX,2 h 

r z max 

dz d ^ 


Jo 

dz 


wv z)w * {zf 


(23) 


with Wf(z) defined as 


p Mm,ax J 

W?(z) = J dM—(M,Z)b(M,z) Sl (M,z) (24) 
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Similarly, the X-ray-galaxy cross-power spectrum is given 
by 

eft# = ^ f dz % f dM m Sab{M ’ z)si{M ^ z) 

eft# = J dz lN p ™{ k = ^- z y)w&(*)wT(*) ( 25 ) 

To calculate < Ssxrb > we use the SXRB 
count s expected to be o bserved by the eROSITA mis¬ 
sion d Merloni et alj 120121 ) assuming a conversion factor 
1 count/sec ~ 10~ n ergsec _1 cm _2 , which is approximately 
the conversion factor for the ROSAT all sky survey. Note 
that using a constant conversion factor underestimates the 
power at low redshifts (z<l for massive galaxies) and over¬ 
estimates the power at high redshifts (z>l). Therefore, this 
simplified approach gives only an order of magnitude es¬ 
timate of the X-ray power spectrum. However, this does 
not affect the estimate of the uncertainty on the model pa¬ 
rameters using X-ray-galaxy cross-power spectrum as the 
main contribution to the constraints comes from the mas¬ 
sive galaxies near z~l. In addition, these constraints satu¬ 
rate fast once the information from two or more appropriate 
mass-redshift bins is combined together as shown in section- 

1641 

In Fig-© we show the X-ray-galaxy cross power spec¬ 
trum for a few mass and redshift bins. For the flat den¬ 
sity and temperature profiles, the shape of the X-ray-galaxy 
cross-correlation power spectrum is analogous to the shape 
of the SZ-galaxy cross power spectrum and has the same 
value of Imax- The contribution of the 2-halo term is neg¬ 
ligible compared to the 1-halo term beyond l = 10 4 for all 
mass and redshift bin combinations considered here. 

In Figure-© we show the X-ray-galaxy cross power 
spectrum for Mh : [10 12 ' 8 ,10 13 ' 0 ](/i -1 Mq) and z : [0.4, 0.5] 
for the density profiles described in section- 13741 Note that 
for all these profiles we have assumed the CGM to be at the 
virial temperature and we truncate the signal at the virial 
radius. The X-ray-galaxy cross-correlation signal increases 
rapidly compared to the SZ-galaxy cross-correlation signal 
with the steepening of the density profile. This is due the fact 
that the X-ray emission is proportional to n 2 and is more 
sensitive to the gas density profile compared to the SZ-effect 
which is proportional to n e . The difference between the pro¬ 
files is now significant even at smaller /-values (~ 3000). 


4.2 Detectability of the X-ray-galaxy 
cross-correlation signal 

As mentioned earlier in section- 13.5.11 the cumulative signal- 
to-noise ratio provides an efficient way of estimating the de¬ 
tectability of a signal for a given survey. The cumulative 
SNR of the X-ray-galaxy cross power spectrum in ab th bin 
is 

(jy) — ^ll'GT,( a b){MlV,(a,b)) 1 Ci'\ a b)5lV (26) 


where the covariance matrix M^ ab ^ for the X-ray-galaxy 
cross-correlation is given by 


Mi 


IV ,(ab,a'b') 


fsky [21 + 1)A/ 


I /yxx /Shh . fhxh Aixh 

X Wz ^l : (ab,a'b')-r'~'l,(ab) , ~'l,(a'b') I 

( 27 ) 



l 


Figure 6. The X-ray-galaxy cross power spectra for : 

[10 12 8 , 1O 13,o ](/i _ 1 M 0 ) and 2 : [0.4, 0.5] for the density profile- 
(a) (red solid line), profile-(b) (green dashed line) and profile-(c) 
(blue dot-dashed line). In the inset we show the corresponding 
density profiles. 


with Ci = Cl + . Assuming that the noise in X-rays is 

dominated by the shot noise, the noise in the X-ray power 
spectrum is 


1 „ „ + ^FWHM \ 

N bg P V 8ln2 ) 


(28) 


where Nb g is the total number of the soft X-ray photons 
collected/steradian by the X-ray telescope and #fwhm is 
the full width half maximum of the beam. 

In practice some of the X-ray background is produced 
by the X-ray emission of AGNs, X-ray binaries, SNe rem¬ 
nants amongst which AGNs are dominant. AGNs are clus¬ 
tered with galaxies and would thus contribute to the X-ray- 
galaxy correlation. This could however be mitigated by first 
removing the fraction of X-ray AGNs which are above the 
detection threshold of the X-ray survey. In addition, AGNs 
have a harder X-ray spectrum than the diffuse circumgalac- 
tic X-ray gas. The photon energy dependence of the cross¬ 
correlation signal can thus also be used to separate the con¬ 
tribution from AGNs. In addition, the signal from AGNs 
would be produced by a 1-source term and a clustering term 
with an angular dependence determined by the point spread 
function (PSF) of the X-ray instrument and the correlation 
function of the AGNs. This specific angular dependence can 
be used to disentangle the contribution from the AGNs and 
from the circumgalactic gas. While a detailed analysis which 
incorporates these mitigating techniques is beyond the scope 
of this paper, we make the optimistic assumption that the 
X-ray noise is not correlated with the distribution of galax¬ 
ies, i.e. Ni h = 0 in Eg 1271 


4.3 X-ray survey 

Since we are interested in the X-ray signal from the galaxies, 
we need an X-ray survey with a small beam size (high res¬ 
olution) and large sky coverage. We consider the eROSITA 
survey for this purpose. The eROSIT A is a future mission 
expected to be launced in 2016 (see iMerloni et al.l (120121 ) 
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Figure 7. Contours of the cumulative signal-to-noise (SNR) for 
the measurement of the cross-power spectrum of the X-ray emis¬ 
sion from the CGM with the distribution of galaxies for the com¬ 
bined eROSITA-DES-like survey. The red (solid), green (dashed), 
blue (dotted) and brown (dot-dashed) lines represent the SNR 6, 
5, 3 and 2 respectively. Note that the upper redshift limit plotted 
here is more than the highest redshift probed by DES. 


for the details of this mission). The specifications of this 
mission are given in Table- The total background ex¬ 
pected in the soft band of the eROSITA is ~ 2 x ICC 3 
counts sec -1 arcmin -2 (I Merloni et al.f2012h . To calculate the 
X-ray-galaxy cross-power spectrum, we consider the combi¬ 
nation of the eROSITA-DES and eROSITA-LSST (Large 
Synoptic Survey Telescope) surveys. 

In Fig-0 we show the contours of the cumulative SNR 
for the X-ray-galaxy cross power spectrum for the eROSITA- 
DES combination. Similar to the SZ-galaxy cross power 
spectrum, the galaxies are binned in mass and redshift with 
A log (M) = 0.2 and A z = 0.1 respectively. Again the 
high mass galaxies have larger SNR due to their larger gas 
reservoir compared to the low mass galaxies and the SNR. 
increases with increasing redshift, becomes maximum and 
then decreases with further increase in redshift. The differ¬ 
ence between the X-ray-galaxy and SZ-galaxy cross power 
spectra is that the X-ray cross-power spectrum peaks at 
relatively smaller redshift compared to the SZ cross-power 
spectrum due to the fact that the observed X-ray surface 
brightness decreases rapidly with increasing redshift. Also 
the contribution from the low mass galaxies at high redshifts 
is more than that of the SZ-galaxy cross power spectrum. 
This is essentially due to the much better resolution of the 
eROSITA (~ 20”) compared to the resolution of the SPT 
(~ 1'). The SNR peaks for the high mass and intermedi¬ 
ate redshift galaxies (z ~ 1). For the most massive galaxies 
(M h = [10 12 - 8 , 1O 13 O ]/i _ 1 M 0 ) at redshift z ~ 1, the SNR is 
~ 7. For the low mass galaxies (< 10 12 h _1 Mg), the SNR is 
low (< 2) at all redshifts. Massive galaxies have significant 
SNR (~ 6) even at redshifts < 0.5, hence the X-ray-galaxy 
cross-power spectrum can be used to study the CGM of 
these systems. 



Figure 8. The cross power spectra of the SZ effect in RJ limit 
and soft X-ray emission from CGM. Here green (dashed) lines, 
blue (dot-dashed) lines and red (solid) lines represent the 1-halo 
term, 2-halo term and the total signal respectively. 


Table 2. Fiducial values and priors on the parameters 


Parameter 

Fiducial value 

Prior-1 

Prior-2 


0.831 

0.013 

0.013 

9m 

0.3156 

0.0091 

0.0091 

n s 

0.9645 

0.0049 

0.0049 

/gas 

0.11 

- 

- 

/Temp 

1.0 

- 

0.25 

Q!gas 

0.0 

- 



The sky coverage of the LSST is 4 times the sky cover¬ 
age of the DES. Since the cumulative SNR oc \Jf s k y , SNR 
for the eROSITA-LSST combination will be twice the SNR 
of the eROSITA-DES combination. Therefore, for the as¬ 
sumed CGM properties, the eROSITA-LSST survey will be 
able to detect the X-ray-galaxy cross-correlational signal 
from the galaxies with a peak SNR ~ 14. Note that the 
estimates of the SNR for both the SZ-galaxy and X-ray- 
galaxy power spectra depend on the size of the mass and 
redshift bins. Therefore, increasing or decreasing the size of 
mass and/or redshift bin also changes the detectability of 
these signals accordingly. 


5 X-RAY-SZ CROSS-CORRELATION 

The SZ effect and X-ray emission have different dependences 
on the gas density and temperature. The SZ effect is pro¬ 
portional to n e T e whereas the X-ray emission scales approx¬ 
imately as n 2 %/Te. Combining the two can improve the con¬ 
straints on the gas physics parameters. The X-ray-SZ cross 
power spectrum is given by 

(jxy _ ^jxy,\h ^jxy,2h 


(29) 


10 0 counts/sec/arcmirf 
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where 




(jxy,lh 

= g(x») [ 

Jo 

. dV f Mma * 

d *ji h 

' lvl min 





xyi(M,z)\ 

Qxy,2h 

= g(xv) [ 

Jo 

dV 

dz — Pm(k = 
dz 

^y^(z) 


xW^z)} (30) 


In Figure-© we show the X-ray-SZ cross power spec¬ 
trum for the flat density profile and the mass and redshift 
range specified in section- l3Tl in dimensionless as well as 
10 _6 counts sec^arcmuG 2 units. The cross-power spectrum 
peaks at l « 15000 for galactic haloes where as for galaxy 
cluste rs it peaks at l ~ 3000 (see Figure-1 of lHurier et al.l 
(l2014h ~). This difference is mainly because the galaxies are 
smaller objects than the clusters and therefore, the X-ray-SZ 
cross power spectrum for galaxies peaks at smal l er ang ular 
scales or larger 1-values. Recently, iHurier et al.l (l2015l f de¬ 
tected the total X-ray-SZ cross-power spectrum at 28cr level 
with ROSAT and Planck all sky surveys. For the Z-range of 
interest for the CGM, we show the SPT-eROSITA combina¬ 
tion. However, due to the weak signal compared to the noise 
for SPT-eROSITA surveys (SNR~ 0.65), it is not possible to 
detect the X-ray-SZ cross-correlation signal from the CGM. 


6 FORECASTING OF CGM CONSTRAINTS 

We now use the Fisher matrix formalism to forecast the 
expected constraints for different survey combinations. Since 
the cosmological parameters are well constrained by Planck, 
our main focus is to constrain the astrophysical parameters 
related to the gas physics. The parameters considered for 
this work are 

{[(T8 , Hm, n s ], [/gas,/Temp, Ogas]} , (31) 

where the parameters in the first bracket are cosmologi¬ 
cal parameters and in the second bracket are astrophys¬ 
ical parameters which depend on baryonic physics. Note 
that we have assumed the flat-A-cold dark matter cosmol¬ 
ogy. The fiducial values of these parameters are given in 
Table-®. Here / gaB is the redshift independent gas frac¬ 
tion, /Temp = tjt— , i.e., the ratio of the temperature of the 
gas to the virial temperature of the gas in the halo and a gas 
represents any possible evolution of the gas defined through 
fsas(z) = fga.s[E(z)] as * s . Our fiducial model assumes no red¬ 
shift evolution of the gas fraction i.e. a gaB =0. 

Given a fiducial model, the Fisher matrix can be written 
as 

Fij = (32) 

dpi dpj 

where M U / is the covariance matrix given by equation- Gi). 
To calculate the uncertainty on the parameters we have con¬ 
sidered following two prior cases: 

Prior-1 : Priors on cosmological parameters only. 

Prior-2 : Priors on cosmological parameters + 25% prior on 

/Temp ■ 

We have included only those galaxies for which the gas 
cooling time is more than the halo destruction time ensuring 
that the CGM temperature is close to the virial temperature 


Table 3. Error on parameters for different scenarios (see Table® 
for SPT-DES combination 


Parameter 

pi 

P2 

PI (fixed ct gaa ) 

P2 (fixed o gaB ) 


0.049 

0.037 

0.042 

0.025 


0.369 

0.207 

0.369 

0.207 

C^gas 

0.519 

0.519 

- 

- 


Table 4. Error on parameters for different scenarios (see Table® 
for eROSITA-DES combination 


Parameter 

pi 

P2 

PI (fixed a gaB ) 

P2 (fixed o gaB ) 


0.20 

0.025 

0.036 

0.015 

/Temp 

2.65 

0.25 

0.649 

0.233 

£*gas 

1.219 

0.318 


" 


of the halo. Therefore, we assume a small uncertainty in 
/Temp in Prior-2 . 

6.1 Results 

In Table-® and @ we show the forecasted uncertainty on 
the parameters for the SPT-DES and eROSITA-DES sur¬ 
veys respectively. Here we have combined the Fisher from 
three different combinations of the mass and redshift bins 
M1Z1, M2Z1 and M2Z2 where Ml=[10 118 ,1O 12 O ]Zi _ 1 M 0 , 
M2=[10 12 ' 8 , lO 13 '°]h _1 M 0 , Zl=[0.4, 0.5] and Z2=[1.0,1.1], 
We show the constraints on the astrophysical parameters 
only as the cosmological parameters are already well con¬ 
strained by Planck. There is a strong degeneracy between 
the gas physics parameters if we consider only one mass and 
redshift bin. However, when the information from two or 
more bins are added together, we can break this degener¬ 
acy and obtain strong constraints on these parameters. In 
Table-© we show the constraints on the parameters from 
the SZ-galaxy cross-correlation signal for the SPT-DES sur¬ 
vey. Combining the Fisher matrix from M1Z1, M2Z1 and 
M2Z2 can constrain / gas to ~ 44% and /Temp to ~ 37% 
around their fiducial values, even without any prior knowl¬ 
edge on astrophysical parameters. For Prior-1 , the constraint 
on a gaB is Aa gas ~ 0.5. Including additional 25% prior on 
gas temperature does not improve the constraint on ct gaB 
whereas the constraint on / gas (/Temp ) improves consider¬ 
ably to 34% (21%). 

In the absence of any redshift evolution of the gas frac¬ 
tion, the constraint on /Temp does not improve whereas the 
constraint on / gaB improves to ~ 38% for Prior-1 and to 
- 23% for Prior-2 . 

We use only above three mass-redshift bins to forecast 
the constraints on gas physics parameters as the addition of 
more bins does not improve these constraints much. For ex¬ 
ample, for Prior-1 , the addition of mass-redshift bin M2Z3, 
where Z3=[0.8, 0.9], improves the constraint on / gaB and a gaB 
from 0.049 and 0.5 to 0.047 and 0.46 respectively whereas the 
change in the constraint on /Temp is <1%. Also the change 
in these constraints in other prior cases is negligible. Further 
addition of mass-redshift bins in the Fisher matrix analysis 
does not improve these constraints. Therefore, we use only 
the bins M1Z1, M2Z1 and M2Z2 for the purpose of our work. 

In Table-© we show the constraints on the parame- 
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Figure 9. The 68% CL for the gas physics parameters / gas , /Temp and o:„ as from the SPT-DES (SZ-galaxy cross power spectrum) and 
eROSITA-DES (X-ray-galaxy cross power spectrum) surveys in the upper and lower panels respectively. All the plots are for Prior-2 case 
and in the right panels, we have fixed a ga s ■ The brown (dot-dashed), green (dashed), blue (dotted) and red (solid) lines represent the 
ellipses for the mass and redshift bins M1Z1, M2Z1, M2Z2 and the combined ellipse respectively. 


ters from the X-ray-galaxy cross-correlation signal for the 
eROSITA-DES survey. Now for Prior-1 and in the presence 
of unknown redshift evolution of the gas fraction, astrophys- 
ical parameters are poorly constrained by this survey. This is 
mainly due to the large noise contamination from the X-ray 
background. The addition of a 25% prior on /Temp signifi¬ 
cantly improves the constraints on the parameters. The un¬ 
certainty on / g as reduces to 23% and a gas can be constrained 
to Aaigas r ^ J 0.3. 

In the absence of any redshift evolution of the gas frac¬ 
tion, / g as (/Temp ) can be constrained to ~ 33%(65%) for 
Prior-1 and the constraint becomes ~ 14%(23%) for Prior-2 

In Figure- ([I)]) we show the constraints from the SPT- 
DES (SZ-galaxy cross power spectrum) and eROSITA-DES 
(X-ray-galaxy cross power spectrum) surveys in the upper 
and lower panels respectively. All the plots are for Prior-2 
and in the right panels, we have fixed a gas . In the upper 
left panel we show the 68% confidence limit (CL) ellipse for 
/gas and agas for the SZ-galaxy cross power spectrum. The 
individual ellipses for M1Z1, M2Z1 and M2Z2 are large and 
there is a large uncertainty on these parameters. But com¬ 


bining them together results in A/ gas « 0.037 i.e. we can 
constrain / gas to ~ 34%. This is because the degeneracy 
of /gas with a gas is broken when we add information from 
the galaxies in similar mass bins but in different redshift 
bins. The X-ray-galaxy cross power spectrum also has simi¬ 
lar contours for / gas v/s a gas (lower left panel of Figure- ©). 
In this case, the Fisher matrix analysis gives a constraint 
A/gas ~ 0.025 on gas fraction, i.e. 23% of its fiducial value. 
In both the cases, M1Z1 bin has relatively large uncertainty 
and the final uncertainty ellipse is essentially determined by 
the other two bins. The amount of gas present in the M1Z1 
bin is roughly an order of magnitude smaller than that of 
the other bins which results in a smaller signal and large 
uncertainty on the parameters. 

In the upper right panel of Figure- © we show the con¬ 
straints on the SPT-DES survey for Prior-2 case and the 
redshift-independent gas fraction. Here adding the informa¬ 
tion from different bins does not improve the constraints on 
the parameters as we already have a strong prior on /Temp 
(~ 25% of its fiducial value). However, even with a strong 
prior on the CGM temperature, there is a large uncertainty 
in the gas fraction for M1Z1 bin which is again due the small 











12 Singh, Majumdar, Nath, Refregier and Silk 



Figure 10. The 68% CL contours for / gaB and 7 gas from the 
SPT-DES (SZ-galaxy cross power spectrum) and eROSITA-DES 
(X-ray-galaxy cross power spectrum) surveys. The solid thin red 
(thick green) ellipse is for the fiducial value 7 gas =0 for the SZ (X- 
ray)-galaxy cross power spectrum and the dashed thin red (thick 
green) ellipse is for the fiducial value 7 gas =2 for SZ (X-ray)- 
galaxy cross power spectrum. All the plots are for Prior-2 case 
with fixed a gas . 


signal in this bin. In the lower right panel on Figure-® we 
show the 68% CL ellipses for / gas -/Temp from eROSITA- 
DES survey. This survey can constrain the gas fraction to 
A/ gas ss 0.015 in case of a redshift-independent / gas . 

Note that, in the calculation of the constraints on 
gas physics parameters, we ha ve used strong prior s on t he 
cosmological parameters (from [Planck results X I ill ( 20151 )). 
However, if we remove the prior on the cosmological param¬ 
eters, it increases the uncertainty in the determination of 
the gas physics parameters. For example, if we remove the 
prior on as, the uncertainty on the gas physics parameters 
becomes larger than 100% for SZ-galaxy as well as X-ray- 
galaxy cross-correlation in Prior-1 case. But the inclusion 
temperature prior (Prior-2 case) improves the constraint on 
/ gas (A/ gas = 0.041 (0.026) for SZ-galaxy (X-ray-galaxy) 
cross-correlation) and a gaB (Aa gaa = 0.52 (0.38) for SZ- 
galaxy (X-ray-galaxy) cross-correlation) significantly even in 
the absence of any prior on as- These constraints are simi¬ 
lar to the constraints on these parameters in the presence of 
Planck prior on as- Therefore, the prior knowledge of one of 
the gas physics parameter helps in breaking the degeneracy 
between the cosmological and gas physics parameters. 


7 PROBING THE ENERGETICS OF THE CGM 

Currently, there is a large uncertainty in the knowledge of 
the distribution of the CGM. Simulations suggest that the 
extent and distribution of the CGM also depends on the 
feedbac k processes taking p l ace in the central p art of the 
galaxy (iLe Brun et alJ 120151 ; ISuresh et al.l 120151) . Without 
any feedback, the temperature at outer radii falls rapidly, 
whereas the winds and AGN feedback tend to make the 
profile flatter. The effect of the feedback on the density pro¬ 


file is however weaker than on the temperature profiles. In 
other words, the pressure profile is likely to be flatter at the 
outer radii in the presence of the feedback processes than in 
the cases without any feedback. 

In this section, we investigate whether the SZ/X-ray- 
galaxy cross power spectrum can distinguish between dif¬ 
ferent density profiles of the CGM and hence the processes 
giving rise to these profiles. In order to constrain the density 
profile of the CGM, assuming the gas to be at the virial tem¬ 
perature, we now include 7 gas in the Fisher matrix analysis, 
where 7 gas is defined by 

Pgas = Poll + ir/RsV^T 1 (33) 

where po is the normalization such that the mass within the 
virial radius of the galaxy remains the same and the fiducial 
value of 7 gas is 0 for the flat density profile. For the SPT-DES 
survey, the uncertainty on the 7 gas is large and it can at best 
be constrained to 7 gas < 3.16 in the absence of any redshift 
evolution of / gas and Prior-2 case. Even when we include the 
redshift evolution of gas fraction and remove the prior on gas 
temperature, the uncertainty on 7 gas degrades only slightly 
to 7 gas < 3.35. This shows that the SZ-galaxy cross power 
spectrum is less sensitive to the density profile of the CGM 
as compared to other astrophysical parameters within the 
resolution of the SPT (Z ~ 10 4 ). Varying the density profile 
affects the SZ-galaxy cross power spectrum only at large Z- 
values and hence this situation can only be improved by a 
higher resolution CMB survey in the future. 

On the other hand, the X-ray emission which is pro¬ 
portional to the square of the density of the CGM, is much 
more sensitive to its density profile even at small Z-values. 
Also the resolution of the eROSITA (Z ~ 30000) is much 
better than the resolution of the SPT. As a result, 7 gaB can 
be constrained to 7 gas < 0.6 by the X-ray-galaxy cross power 
spectrum. This constraint remains almost invariant even if 
we fix Q gas in the Prior-2 case. 

In the case of a steeper density profile i.e. a larger value 
of 7gas , both the SZ-galaxy and X-ray-galaxy cross power 
spectrum increase, specially at large Z-values (see Figure- © 
and ©). This results in an increased signal-to-noise ratio, 
as well as an improvement in the constraints on the param¬ 
eters. For example, if we take 7 gas =2, the constraint on 
7gas improves to A 7 gas < 1 from the SZ-galaxy cross power 
spectrum whereas the slope of the density profile can be con¬ 
strained with an accuracy better than 5% (i.e. A 7 gaB < 0.1) 
from the X-ray-galaxy cross power spectrum. 

In Figure- (El, we show the 68% CL ellipses for / gaB 
-7gas from the SPT-DES and the eROSITA-DES surveys. 
Again, due to the more sensitive dependence of the X-ray 
emission on the CGM density profile as compared to the 
SZ-effect, the uncertainty ellipses of the X-ray-galaxy cross¬ 
correlation are smaller and hence can constrain 7 gas better 
than the SZ-galaxy cross-correlation. 

This result has implications of being able to constrain 
the di stribution and evolution of the CGM. R ecent simula¬ 
tions (Ihe Brun et alJ 120151 : ISuresh et al.l 120151 ) have shown 
that the feedback processes (from star formation and AGNs) 
can affect the density and temperature profiles of the CGM. 
These simulations match some of the observed properties 
of the CGM and galaxies. However, they are quite sensi¬ 
tive to the feedback mechanism used and give a variety of 
the CGM density and temperature profiles depending on 
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the feedback process. The density profiles of the CGM in 
these simulations can be reasonably fit by 1.2 Sj 7 gaa < 2.5, 
within the virial radius of the galaxy for various feedback 
processes. For example, the No AGN, fixed-v hot w inds and 
fully enriched winds models of ISuresh et alj (l2015l l can be 
fit with 7 gaa =2, 1.6 and 2.2 respectively (excluding the cen¬ 
tral part). Also the pressur e profile for the massive galaxies 
from Le^run_e£jd) (2015) (see the first panel of Figure- 
(3) of lLe Brun et all (j 20151 )) can be fit with a 7 gaa ~ 1.25, 
assuming the gas to be at the virial temperature. There¬ 
fore the X-ray-galaxy and SZ-galaxy cross power spectrum 
have the potential of discriminating between the evolution¬ 
ary processes for the CGM, at greater than 3a, if A 7 gaa can 
be constrained within ~ 0.5. 


8 CONCLUSIONS 


(SZ)-galaxy cross power spectrum for the flat density pro¬ 
file. These constrains are sensitive to the fiducial value of 
7 gas and improve for a steeper density profile of the CGM. 
For 7 gas =2, it can be constrained to A 7 gas < 0.1 (< 1.0) 
by the X-ray (SZ)-galaxy cross power spectrum power spec¬ 
trum. 

(vi) In all our calculations, we have assumed / gas =0.11. 
If instead, we take a low value of the gas fraction in the 
CGM, e.g. / gaa « 0.05, the SZ cross-power spectrum and its 
detectability goes down roughly by a factor 2 as the SZ signal 
is proportional to the amount of gas. So with / gaa « 0.05, 
the peak SNR ~ 4-5 for the SPT-DES survey. However, since 
the X-ray signal is proportional to , the X-ray cross power 
spectrum and its detectability goes down by a factor of 4. 
As a result, the SNR goes below 2 for the eROSITA-DES 
survey whereas this signal can still be detected at SNR ~ 3 
with the eROSITA-LSST combination. 


We have studied the cross-correlation power spectra for the 
SZ-galaxy distribution, X-ray-galaxy distribution and X- 
ray-SZ effect for the hot gas in the galactic haloes. Our main 
conclusions are as follows: 

(i) We predict that the SZ-galaxy cross power spectrum 
is significant at small scales (l > 3000) and can be detected 
at SNR ~ 9 by combining the SPT and DES surveys for 
the massive galaxies at intermediate redshifts. The shape 
of the SZ-galaxy cross-power spectrum is sensitive to the 
underlying distribution of the gas at large Lvalues (l « 10 4 ). 

(ii) For the X-ray-galaxy cross power spectrum, we have 
considered the combination of the eROSITA-DES and 
eROSITA-LSST surveys and these surveys can detect the 
signal at SNR~ 6 and 12 respectively for the high mass and 
intermediate redshift galaxies. For the flat density profile, 
the shape of the X-ray-galaxy cross power spectrum is sim¬ 
ilar to the shape of the SZ-galaxy cross power spectrum. 
However, the X-ray emission (oc n a ) is more sensitive to the 
density profile than the SZ-effect (oc n e ). As a result, the 
X-ray-galaxy cross power spectrum vaires significantly with 
the steepening of the density profile even at l « 3000. 

(iii) The possibility of detecting the X-ray-SZ cross power 
spectrum from the CGM is low (SNR < 1) for the SPT- 
eROSITA combination. This is due to the combined effect 
of the high noise in X-rays and low resolution of the SPT 
survey. 

(iv) Finally, we do a Fisher matrix analysis for these sur¬ 
veys to forecast the constraints that can be derived on the 
amount of gas in the CGM. After marginalizing over the cos¬ 
mological parameters with Planck priors and combining the 
Fisher matrix analysis for three different mass and redshift 
bin combinations, the SPT-DES survey can constrain / gaa to 
~ 34% in the presence, and to ~ 23% in the absence, of any 
possible redshift evolution of the gas fraction. For the same 
set of mass and redshift bins, the eROSITA-DES survey can 
constrain / gas to ~ 23% and ~ 14% in the presence and ab¬ 
sence of redshift evolution of gas fraction respectively. Note 
that we neglect the correlation between the galaxies and the 
AGNs in the calculation of the uncertainties in the X-ray- 
galaxy cross-power spectrum. 

(v) Including the slope of the density profile 7 gas (defined 
in Ea 1331) in the Fisher matrix analysis, degrades the con¬ 
straints on other astrophysical parameters whereas 7 gas it¬ 
self can be constrained to 7 gas < 0.6 (< 3.4) by the X-ray 


Presently, the amount of CGM in galactic haloes, its dis¬ 
tribution, energetics and other properties are not well deter¬ 
mined. Therefore, the detection and study of the SZ-galaxy 
and X-ray-galaxy cross power spectrum can provide power¬ 
ful constraints on the nature of the CGM and open up the 
possibility of differentiating between various feedback mod¬ 
els which affect the evolution of the CGM. 
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